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Transport through an Anderson junction (two macroscopic electrodes coupled to an Anderson 
impurity) is dominated by a Kondo peak in the spectral function at zero temperature. The exact 
single-particle Kohn-Sham potential of density functional theory reproduces the linear transport 
exactly, despite the lack of a Kondo peak in its spectral function. Using Bethe ansatz techniques, 
we calculate this potential exactly for all coupling strengths, including the cross-over from mean- 
field behavior to charge quantization caused by the derivative discontinuity. A simple and accurate 
interpolation formula is also given. 



It is a truth universally acknowledged that many-body 
effects in strongly correlated systems are not reproduced 
by mean-field theory. Although Kohn-Sham (KS) density 
functional theory (DFT) is formally exact, it is a non- 
interacting theory yielding only the ground-state energy 
and density of a system. No other information about the 
correlated many-body wavefunction is available. Dynam- 
ical properties, such as excitations and response func- 
tions, are also not predicted by ground-state DFT, even 
with the exact functional [1]. The hope is that, for 
weakly-correlated systems in which ground-state DFT 
approximations perform well for total energies, geome- 
tries, etc., the errors in such calculations are small. Noth- 
ing in the theorems of DFT guarantees that a ground- 
state KS calculation can describe transport correctly [2]. 

Consider transport through an Anderson junction 
[3, 4], composed of two macroscopic leads coupled to an 
Anderson impurity. As an integrable system, the An- 
derson model is a paradigm of many-body physics. It 
is also an accurate model of the low-energy spectrum of 
molecular radical-based junctions [5]. In Fig. 1, we show 
the zero-temperature transmission through an Anderson 
junction as a function of the energy e of the resonant level 
using Bethe ansatz (BA), exact Kohn-Sham (KS) DFT, 
and Hartree-Fock (HF). In the figure, fi is the chemical 
potential (Fermi energy) of the metal electrodes. Re- 
markably, the KS-DFT treatment of this problem repro- 
duces the transmission exactly, apparently describing the 
non-pcrturbative Kondo effect whose spectral peak is the 
source of the perfect transmission when e < fi < e + U . 

The inability to describe sharp steps in transmission is 
a well-understood failure of standard density functional 
approximations. In the limit of weak coupling to the 
leads, the system is a prototype example where the ef- 
fects of the infamous derivative discontinuity is seen [6] . 
For such a system, the exchange-correlation (XC) energy 
of the molecule is strictly linear between integer values, 
and so the XC potential, its functional derivative, jumps 
discontinuously at such values [6]. This effect has been 
implicated in many well-known failures of DFT approx- 
imations such as strongly-correlated systems [7], charge- 
transfer excitations [8], and over-estimation of the cur- 
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1. (Color online) Zero-temperature transmission of an 
Anderson junction as a function of e using Bethe ansatz 
(exact), Kohn-Sham DFT with the exact functional (KS), 
and (spin-restricted) Hartree-Fock (HF). As U increases, HF 
misses the sharp structure, but the KS transport is always 
exact! 

rent in organic junctions [9]. 

In this letter, we (a) solve the Anderson junction us- 
ing BA and invert the KS equations to derive the ex- 
act KS potential, (b) show that the transport is recov- 
ered exactly, but only for zero temperature and weak 
bias, and (c) parametrize the XC potential, providing a 
unique interpolation formula that works for all correla- 
tion strengths. The exact and KS-DFT transport are 
intimately related via the Friedel sum-rule, a statement 
equating the occupancy and the transmission phase at 
the Fermi energy, as has been discussed for nearly iso- 
lated resonances [10, 11] or single-mode leads [11]. Since 
our system fits neither of these categories, it reopens the 
question of just when an accurate ground-state KS cal- 
culation will yield accurate transport. 

Using BA [12, 13], we calculate the exact KS potential 
for the Anderson model, and investigate how the deriva- 
tive discontinuity develops in the limit of weak impurity- 
lead coupling. From this solution, we extract the exact 
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asymptotic scaling form of the derivative discontinuity 
and establish an interpolation formula for the KS po- 
tential which is accurate for all coupling strengths, com- 
putationally simple, and illustrates the crossover from 
mean-field behavior at strong impurity-lead coupling to 
charge quantization at weak coupling. 

Consider a junction composed of a nanoscale central 
region (C) connected to two macroscopic electrodes, la- 
beled left (L) and right (R) . The Hamiltonian of the sys- 
tem is % = He + Ht + Hr + Hl, where "Hr/l describe 
Fermi gases in the electrodes and Ht describes tunnel- 
ing between the central region and the electrodes. The 
central region has the Hamiltonian [3]: 



(1) 



where n a = d\d a is the number operator for spin a and 
the charging energy U is given by the Coulomb integral. 
The exact Green's function of an Anderson junction may 
be found using Dyson's equation in an orthonormal basis: 



G(E) = [E-e- E C (E) - E T {E)Y 



(2) 



where Sc is the Coulomb self-energy and St 
is the tunneling self-energy [14]. We consider transport 
in the broad-band limit where ~E T (E) = —iT a /2 is pure 
imaginary and independent of energy, and define the 
mean tunneling-width T = (T L + T R )/2. 

At zero temperature, the linear-response transmission 
function of an Anderson junction is given by [4, 13] 



T(£)| B=M = r L r R |G( M )| 2 = 



4PT 



Lt^R 



(r L + r R )' 



sin 



9[p) 



where 0((J.) is the sum of transmission eigenphases eval- 
uated at the Fermi energy /i = /i^ = /i^j. The total 
number of electrons on the central impurity is related to 
9{ji) at zero temperature by the Friedel sum-rule [15-17] 



(n c ) = 0(/i)Ar. 
From the Bethe ansatz, one finds [12] 

i f°° duj e^l/ 2 
(no(M)>=l- ' 



e ^(»(fc)-Q)A(fc)dfc 



(4) 



(5) 



with Q defined by the condition 

U - 2fJ, _ i f°° doj e -M/2-i"Q 

VWr ~ 7^ J-oo u + iv g(-)(w) ^f- 

making use of the following functions 



2tt 



r(l/2 + io;/27r) 



n 



2ire 



(6) 



(7) 



g(k) = (k + ^-U/2) 2 /(2UT), and A(fc) = (l/vr)r/(r 2 - 
(k+/j,) 2 ), where T(x) is the Gamma function and rj = 0^ 
(nc) is plotted as a function of e in Fig. 2. 
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FIG. 2. (Color online) The occupancy of the Anderson junc- 
tion as a function of fi — e. The BA results (solid line) exhibit 
discrete charge steps when U 3> T, whereas the mean-field 
RHF (dashed line) never does. The approximate occupancies 
calculated self-consistently within the KS-scheme using exc 
(Eq. (15)) (dotted line) are nearly indistinguishable from the 
exact BA results. The exact occupancies were used to gener- 
ate Fig. 1. 

The KS ansatz of DFT employs an effective single- 
particle description, defined to reproduce the ground- 
state density of the interacting system. By the 
Hohenberg-Kohn theorem, this potential is unique if it 
exists (it usually does) [1]. The relationship between po- 
tential and density is fixed only in the full basis-set limit, 
but can be defined for lattice models. The leads of an 
Anderson model are non-interacting and characterized 
by a total charge, which remains constant. Thus the def- 
inition of the KS system in this extreme case is that of 
a non-interacting junction (U = 0) with an on-site po- 
tential chosen to reproduce the on-site occupancy of the 
interacting system. 

The KS Green's function in the central region may be 
written as 



G S {E) = (E-e s + iTY 



(8) 



where e s is the KS potential for an electron on the impu- 
rity, and is written 



e s =e + U(n c )/2 + exc, 



(9) 



where the second term is the Hartrce contribution and 
the last is the correlation potential (there are no ex- 
change contributions) . The Anderson model has no inter- 
nal molecular structure so the KS lead-molecule coupling 
is r [5, 18]. For more complex systems, T need not be 
equal to the KS lead-molecule coupling. In a standard 
DFT calculation, £xc is approximated as a functional of 
the density [1]. The occupancy of the central region is 

^Im[G(E)}<, (10) 



nc = 2 



where the "lesser" Green's function is found using the 
Keldysh relation [14] 



[G(E)] < = 2if(E)T\G(E)\' 



(11) 
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FIG. 3. (Color online) Transmission through an Anderson 
junction at fixed n — 3eV with e = 0, V = 0.5eV and 
U = WeV, so that T K = 2mK [12, 13, 19]. The Doniach- 
Sunjic form [20] of the spectral function [21] is used near the 
Kondo peak, while the non-singular portion is calculated us- 
ing the methods of Refs. [5, 14] with the exact occupancy. 
Logarithmic shifts of the charging resonances are also included 
[19]. More sophisticated numerical methods [22] reproduce 
the qualitative features shown here. 

where at zero temperature, f(E) = 0(/i— i?) and is the 
Heaviside function. Inserting the KS Green's function 
and solving Eq. (10) for e s gives 



e s = /i + rcot(|(n c )) 



(12) 



which defines the exact KS potential within the Anderson 
model. The KS transmission is then 



T S {E) 



FT 



{E-e s y + r 2 



(13) 



Plugging Eq. (12) into Eq. (13), we find that T S (E) is 
identical to T(E), as was shown in Fig. 1. Although this 
identity can be derived using, e.g., local Fermi liquid the- 
ory [4], nonetheless its significance is profound: If (and 
only if) a mean-field theory yields the exact occupation 
will it yield the exact transmission. 

In an Anderson junction, the Friedel sum-rule connects 
the transmission at the Fermi energy to the occupancy 
at zero temperature. As shown in Fig. 3, the full trans- 
mission spectrum of an Anderson junction exhibits three 
peaks: Two Coulomb-blockade peaks of width ~ r cen- 
tered about E w and E « U and a third zero-bias 
Kondo peak of width ~ ^bTk pinned at E = fx. In 
contrast, the KS-DFT transmission spectrum is a single 
Lorentzian of width T peaked at E = [i. As indicated 
in the figure, the KS value is a huge overestimate any- 
where more than several IcbTk away from u, implying 
that the ground-state KS potential does not accurately 
predict transport at temperatures larger than Tr or for 
bias voltages larger than fceTft-/e. 

From Eqs. (3), (4), (12) and (13) it is evident that the 
HF errors in transmission in Fig. 1 stem from errors in 
the occupancies of Fig. 2. When U < T, HF yields accu- 
rate occupancies and transmissions. But for U ^> T, the 
HF occupancies lack the distinct steps present in the ex- 
act solutions, causing corresponding discrepancies in the 
transport. Qualitatively similar errors would be found 
with any local or semilocal approximation for the XC 



potential, because such approximations are smooth func- 
tions of the interaction strength [23]. But the exact KS 
potential of an isolated system, infinitely weakly coupled 
to a reservoir, displays discontinuous jumps at integer 
particle number [6], just as ours does as T/U — > 0. The 
KS potential e s is shown as a function of occupancy in 
Fig. 4a for several values of U/T. The HF potential is 
linear with a slope of U/2. For large but finite U/T, the 
KS potential is not discontinuous but has steps (of width 
~ T/U) corresponding to the plateaus in the occupancy, 
becoming discontinuous in the limit. When U/T is suffi- 
ciently small there is no step in the KS potential and the 
HF approximation is accurate. 

We now show how the step develops as T/U — > 0. In 
Fig. 4a, e s develops a step of height U at (nc) = 1 whose 
width decreases as T/U — > 0. Fig. 4b shows the exact 
derivative ds s /d(nc) in the vicinity of (nc) = 1. The 
horizontal and vertical axes are rescaled to illustrate the 
scaling behavior of the step as T/U — > 0. From the BA 
solution, as U — > oo [24]: 



£xc 



U 
2~ 



1 - (nc) - 



tan 



ir 2 U(l - (n c )) 
8T 



(14) 

whose derivative yields a Lorentzian. In Fig. 4b, we show 
this limit and how it is approached as U grows, but notice 
also that the Lorentzian shape is approximately correct 
for all U down to T. We thus parametrize the XC poten- 
tial for U > r by the approximate form 



U 



1 



(n c ) 



tan 



1 



{nc. 



, (15) 



where a and a are functions of T/U which — > 1 and 
8T/(n 2 U), respectively, in the limit T/U 0, and de- 
termine the amplitude and width of the correlation con- 
tribution as a function of (nc). Varying a between 
and 1 corresponds to "turning on" charge quantization, 
a nonperturbative interaction effect that would not be 
described by any local or semilocal approximation. 

To determine those parameterizations, we first note the 
behavior near (nc) = 0, 2 for T ^ 0. This is not generic, 
but stems from the restricted Hilbert space of the Ander- 
son model. For U/T — > oo, Eqs. (12) and (10) imply that 
de s /d(n c ) - r7r/(2sin 2 [7r(n c )/2]) near (n c ) = 0. How- 
ever, 3 x 10 _5 [l/(nc) 4 ] fits the data most accurately over 
the range 0.1 < U/T < 500 . Precisely the same form is 
applied at (nc) = 2. Then a least-squares fit of Eq. 
(15) to the exact BA results, subtracting the features at 
(nc) =0,2 and limiting the fit range to 0.1 < (nc) < 1.9, 
yields the values of a and a given in Fig. 4c by the points. 

We also found simple fits for these functions, with 



U 



U + 5.68T 1 

r 



(16) 



0.811 



U 



0.390 



IP 



0.168 



U 3 ' 



as shown by the smooth curves in Fig. 4c. The corre- 
sponding derivatives are given by dashed curves in Fig. 
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FIG. 4. (Color online) The exact KS potential e B of the Anderson junction, (a) For U«r, HF is accurate; as U grows, an 
increasingly sharp step develops, due to the derivative discontinuity of the exchange-correlation energy as a function of particle 
number, (b) de s /d{nc) for several large values of U/T. As U/T — > oo, the exact (Lorentzian) asymptotic scaling form is 
recovered, (c) Exact and interpolated a and a functions are shown in the left and right-hand panels, respectively, showing the 
crossover from the meanfield (U/T ~ 0) to the strongly correlated regime (U/T — ¥ oo). The crossover from weakly to strongly 
correlated occurs for U/T « 6. (d) Exact and approximate de s /d{nc) spectra for several moderate U/T values, highlighting 
the accuracy of Eq. (15). The dashed lines use Eq. (16). 



4d. Eqs. (15)— (16) define an interpolation formula for 
the XC potential which yields the exact KS potential in 
both the U/T — > and U/T —> oo limits and is accu- 
rate for all intermediate values of U > T, as shown in 
Fig. 4. To check that our parametrization is sufficiently 
accurate, we performed self-consistent calculations using 
Eqs. (15) and (16), finding transmissions indistinguish- 
able from the exact ones for all values of U/T (see Fig. 2). 
Our interpolation formula shows that the cross-over be- 
tween weak and strong correlation occurs for U ~ 6r. 

Our results should prove useful for the development of 
density functional theory in general, and for its appli- 
cation to transport through molecular junctions. While 
much is known of the consequences of the derivative dis- 
continuity in the extreme limit of weak coupling, our re- 
sults describe the entire range from strong to weak cou- 



pling, i.e., from weak to strong correlation. For transport 
through molecular junctions, our results provide an exact 
limit for which both many-body and DFT approxima- 
tions can be tested. Nor are these just theorists' games 
with toy models. For example, for the archetypal Au- 
[l,4]benzenedithiol single-molecule junction, U > 8T. In 
such a system, any approximate XC functional which fails 
to account for derivative discontinuity effects is unlikely 
to yield accurate results. 

Note added in proof: After submission of our 
manuscript, Refs. 18, 25, 26 (chronological order) re- 
porting related results appeared on arXiv. This work 
was supported by DOE under grant number DE-FG02- 
08ER46496. KB thanks David Langreth for a decade- 
long conversation on the subject, and a lifetime of inspi- 
ration. 
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